function H = lowpassFilter(fftStruct, cutoffFreq)

% H: Filter transfer function
[H r c] = initFilterSize(fftStruct.fftImage);

%Cut off Frequency
Do = cutoffFreq; 

for i = 1 : size(H,1)
    for j = 1 : size(H,2)
        dist = sqrt((i - (r/2))^2 + (j - (c/2))^2);
        H(i,j) = dist <= Do;
    end
end

%[u,v] = meshgrid(1:1:r,1:1:c);
%surf(u,v,H);
%axis([-500,500,-500,500,0,1]);

end